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We propose a Bell inequality for high-dimensional bipartite systems obtained by binning local 
measurement outcomes and show that it is tight. We find a binning method for even d-dimensional 
measurement outcomes for which this Bell inequality is maximally violated by maximally entangled 
states. Furthermore we demonstrate that the Bell inequality is applicable to continuous variable 
systems and yields strong violations for two mode squeezed states. 
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The incompatibility of quantum non-locality with 
local-realistic (LR) theories is one of the most remark- 
able aspects of quantum theory. LR theories impose con- 
straints on the correlations between measurement out- 
comes on two separated systems which are described by 
Bell inequalities (Bis) Q. It was shown that Bis are vi- 
olated by quantum mechanics in the case of entangled 
states. Therefore Bis are of great importance for under- 
standing the conceptual foundations of quantum theory 
and also for investigating quantum entanglement. Since 
the first discussion of quantum non-locality by Einstein- 
Podolski-Rosen (EPR) a great amount of relevant work 
has been done and numerous versions of Bis have been 

proposed 0, fl S, 1 S 1 0, 1, m ei m m d ■ 

For bipartite 2-dimensional systems the CHSH Bell- 
type inequality 0] has the desirable property of only be- 
ing maximally violated for a maximally entangled state. 
The CHSH inequality divides the space of correlations be- 
tween measurement outcomes by defining a hyperplane. 
Since a facet of the polytope defining the region of LR 
correlations lies in this hyperplane the CHSH inequality 
is tight. This means that any violation of LR theories oc- 
curring on this particular facet is indicated by the CHSH 
inequality Tightness is a desirable property since 
only sets of tight Bis can provide necessary and sufficient 
conditions for the detection of pure state entanglement. 
There are still many open questions regarding the gen- 
eralization of Bis to complex quantum systems Q ■ For 
example, Bis for bipartite high-dimensional systems as 
e.g. that proposed by Collins et al. [5( are either not 
maximally violated by maximal entanglement [(| or as in 
the case of Son et al. Q were shown to be non-tight [1] . 

In the case of continuous variable systems there is so 
far no known BI formulated in phase space which is max- 
imally violated by the EPR state - the maximally entan- 
gled state associated with position and momentum Q. 
Although Banaszek and Wodkievicz (BW) showed how 
to demonstrate non-locality in phase space II CM their BI 
is not maximally violated by the EPR state [ll| . Another 
approach using pscudospin operators was shown to yield 
maximal violation for the EPR state [12]. However, find- 



ing measurable local observables to realise this approach 
is challenging. Due to the lack of any known BI pro- 
viding answers to these questions we still have no clear 
understanding of nonlocal properties of high-dimensional 
systems and their relation to quantum entanglement. 

In this paper we present a BI for even d-dimcnsional bi- 
partite quantum systems which, in contrast to previously 
known Bis, fulfills the two desirable properties of being 
tight and being maximally violated by maximally entan- 
gled states. These properties are essential to investigate 
quantum non-locality appropriately and for consistency 
with the 2-dimcnsional case. We call Bis fulfilling these 
properties optimal Bis throughout this paper. Then we 
extend optimal Bis to continuous variable systems and 
demonstrate strong violations for properly chosen local 
measurements. 

Optimal Bell inequality-We begin by briefly introduc- 
ing the generalized formalism for deriving Bis for arbi- 
trary d-dimensional bipartite systems Q. Suppose that 
two parties, Alice and Bob, independently choose one 
of two observables A\ or A^ for Alice, and B\ or Bi 
for Bob. Possible measurement outcomes are denoted 
by k a for A a and l a for B with a, b = 1,2, where 
k a ,h S V = {0, l^...,d— 1}. A general Bell function 
is then written as 



(1) 



d-l 



B= 2J X! e ab(k a Jb)Pab(k a Jb), 



a. 6—1 k a ,/(,— 



where P ab (k a ,lb) is the joint probability for outcomes 
k a and lb, and e a {,(fc a ,ib) are their weighting coefficients 
(here assumed to be real). For local-realistic (LR) sys- 
tems each probabilistic expectation of B is a convex com- 
bination of all possible deterministic values. It can thus 
not exceed the maximal deterministic expectation value 
given by 



23^ x = rnax{ £ e ab (k a ,l b )V 

^a,b=l ' 



(2) 



where C = {(ki,k2,h,l2)\ki,k2,h,h G V] is the set of 
all possible outcome configurations. A quantum state 



violates local realism if its expectation value exceeds the 
bound B™£ x . The flexibility in choosing the coefficients 
£ab{ka, h) allows the derivation of all previously known 
Bis [8[ from Eq. (JTJ), e.g. those proposed by Collins et al. 
Q and by Son et al. Q ■ Moreover, we can construct new 
Bis by properly choosing coefficients e a b{k a , h). Our aim 
is to find optimal Bis which fulfil the following conditions: 

(CI) The BI is tight i.e. it defines a facet of the poly- 
tope separating LR from non-local quantum regions in 
correlation or joint probability space. 

(C2) The BI is maximally violated by a maximally 
entangled state. For each bipartite d-dimensional max- 
imally entangled state there exists a basis \j) with 
j = 0, • • • , d — 1 in which this state reads \ipf ax ) = 

As a general method, one could choose the coefficients 
£ab(k a ,lb) freely and examine whether the resulting BI 
satisfies the conditions (CI) and (C2). Here we instead 
propose a method which restricts this choice and is guar- 
anteed to give tight Bis. We assume that the coefficients 
are products of arbitrary binning functions defined by 
each party as 



(a) 



(Tl) 



(T2) 



(T3) 



( R (k) = 



+1 if outcome k G R, 
— f otherwise, 



(3) 



where R is an arbitrarily chosen subset of all possible 
outcomes, i.e. R C V. The coefficients are then given by 



en = Cri(&i)C,SiGi)> e 12 

£21 = CR 2 ( k 2)Cs 1 (ll), £22 = 



-Cfl 2 ( fc 2)Cs 2 (^), 



(b) 




FIG. 1: (a) The coefficient distributions of Bis (Tl), (T2), 
and (T3) for d = 8 are shown with weighting +1 (white for 
en, £12, £2i and grey for £22) and -1 (grey for en, £12, £21 and 
white for £22) in outcome space, (b) Quantum expectation 
values Bqm of (Tl), (T2), and (T3) are plotted. As d in- 
creases, the expectation value of (Tl) reaches the bound 2\/~2 
(solid line), while that of (T2) approaches 2.31 < 2v2 and 
that of (T3) decreases below the local-realistic upper bound 
2 (dashed line). 



(4) 



where R a and Sb are subsets of the outcomes of A a and 
Bt, respectively From Eq. © we find the LR upper 
bound Bg£* = 2. 

We first show that any BI derived by this method is 
tight. The extremal points of the polytope separating 
LR and non-local quantum mechanical correlations are 
associated with all deterministic configurations C . They 
are described by 4<i 2 dimensional linearly independent 
vectors G ku k 2 ,h,h = ( e fci ® ejj © (e fcl ® e h ) © (e fc2 © 
ei ± ) © (efc 2 © e; 2 ), where is the d-dimcnsional vector 
whose fc-th component is 1 and all other components are 
zero. The interior points of the polytope are given by 
convex combinations of these extremal points and repre- 
sent the region accessible to LR theories. We now only 
consider extremal points associated with configurations 
giving the maximal LR value B™£* and denote their num- 
ber by M. For a polytope defined in Ad 2 dimensions at 
least Ad(d — 1) linearly independent vectors are required 
to define a facet. Therefore, if M > Ad(d — 1) the ex- 
tremal points yielding Bf^ define a facet of the poly- 
tope distinguishing LR from non-local quantum mechan- 
ical correlations We assume the number of elements 
in the sets R\, R2, S±, S2 to be n\, 712, mi, rri2, re- 
spectively, where < ni,n2,mi, mi < d—1. We then 
count the number of configuration giving Bf^ x and find 
Al = d 2 (d 2 —d{ni +mi) +n\[mx +7722) + 712 (mi — 7712)) > 



4d(d— 1). Therefore all Bis obtained by this method are 
tight, i.e. they satisfy condition (CI). Note that any loss 
of elements in binned subsets may cause them to become 
non-tight. 

We now discuss the maximal violation condition (C2) 
by considering three different tight Bis obtained via the 
above method. The corresponding choices of the coeffi- 
cients tab for d = 8 are schematically shown in Fig.[lja): 
(Tl) is the sharp binning type which can be realised if all 
outcomes arc identifiable with perfect measurement reso- 
lution. The elements of the subsets arc given as the even- 
numbers, i.e. R\ = i?2 = Si = 32 = {0,2,4, ...} so that 
the coefficients e ao have an alternating weight +1 or —1 
when an outcome changes by one. (T2) is associated with 
unsharp binning resolution and can be used to model im- 
perfect measurement resolution. The subsets are chosen 
as Ri = R 2 = Si = S 2 = {Vk\k = 0, l(mod 4)} where 
k = 0, l(mod 4) indicates that k is congruent to or 1 
modulo 4. The coefficients e a b alternate between +1 and 
— 1 for every 2 outcomes. In type (T3) the measurement 
results are classified into two divided regions by the mean 
outcome [d/2], where [x] denotes the integer part of x. 
The subsets are chosen as Ri = R2 = Si = S2 = {V/c|0 < 
k < [d/2]}. These three types of binning correspond to 
different capabilities in carrying out measurements on d- 
dimcnsional systems. Their properties will yield useful 



3 



insights for testing Bis in high dimensional and continu- 
ous variable systems. 

We examine quantum violations of (Tl), (T2) and (T3) 
by the maximally entangled state \ip™ ax ) with increasing 
dimension d. The measurements A a and Bb are per- 
formed in the bases \a,k) = (l/Vd) J2lZo uj^ k+a ^ \j) 



and \b,l) = (l/\/5)E 



d-l 



j), obtained by quan- 



tum Fourier transformation and phase shift operations 
on Here uj = exp(2iri/d), and a a and (3b are phase 
factors differentiating the observables of each party A a 
and Bb, respectively. The expectation value of the Bell 
function is then given by 



a 



QM 



2 d-l 



e a b(k,l) 



(5) 



As shown in Fig. [TJb), the expectation values of 
(Tl) for even-dimensions are 2v / 2, and those for odd- 
dimensions tend towards 2^/2 with increasing d. This 
is the upper bound for quantum mechanical correla- 
tions which we show by defining a Bell operator as 
& = E a , b Ek,l^b(k,l)\a,k)(a,k\ ® \b,l)(b,l\. From 

Eq. J!]), B 2 = At d ®t d + [A, A] ® [Q%,Ql] where 

Pa = E*C«.(*)|a,*)(a,*l> = £ { Csr t (0IMKM 
and Id is the d-dimensional identity operator. Since 

l|[A,A]|| < llAAll + ilAAll < 2|jA||!lA!l = 2 and 

likewise for || [Q 2 , Qi] || where || • || indicates the supre- 
mum norm, we finally obtain \\B 2 \\ < 8, or ||23|| < 2y/2. 

We calculate the quantum mechanical expectation 
value of B for BI (Tl) by writing the coefficients as 
en = ei2 = £21 = (-l) k+l and e 22 = -(-l) fe+i . For even 
d we use ESio(-l) fe+i sin[5(fc + I + a a + (3 b )} = 
cos7r(a a + /3b) and find the expectation value 

Bqm = cos7r(ai + /3i) + cos7r(ai + /3 2 ) 

+ cos7r(a2 + Pi) - cos7r(a 2 + /3a). (6) 

This expression also holds approximately for sufficiently 
large odd d. Thus we obtain £>qm = 2\/2, i.e. the max- 
imal quantum upper bound, for a\ — 0, a 2 = 1/2, 
0x = -1/4, and (3 2 = 1/4. Figure QJb) also shows the 
maximal expectation values of (T2) which are smaller 
than 2y2 and approach w 2.31 with increasing d. The 
maximal expectation values of (T3) decrease below the 
local-realistic upper bound 2 with increasing d. 

These results show that (Tl) is an optimal BI for even 
d which satisfies conditions (CI) and (C2). The optimal 
correlation operator is then written as E a b = II a ® lib- 
with the local measurement fl a = Efc=o( — l) fe | a J k)(a, k\. 
Finally, we obtain the optimal BI 



B = E n + E 12 + E 21 - E 22 < 2, 



(7) 



where E ab = (E ab ) = J2kdh l ) k+lp ab(k,l) is the corre- 
lation function. Note that for d = 2 Eq. ([7]) is equivalent 
to the CHSH inequality [|[. We have thus shown that 



the perfect sharp binning of arbitrary even dimensional 
outcomes (Tl) provides an optimal BI, while the other 
binning methods (T2) and (T3) tend to neglect quan- 
tum properties and do not show maximal violation for 
maximally entangled states. 

Continuous variable systems-We extend the optimal 
BI (Tl) to a continuous variable system and calculate 
its violation by a two-mode squeezed state (TMSS). This 
state can, for instance, be realised by non-degenerate op- 
tical parametric amplifiers [l3| in photonic systems. It 
is written as |TMSS) = sechr E^Lo tanh n r\n, n) where 
r > is the squeezing parameter and |n) are the num- 
ber states of each mode. In the infinite squeezing limit 
r — > 00, this becomes the normalized EPR state flO ]. 

When directly following the procedure of the finite 
dimensional case two problems arise: First, we obtain 
the local measurement basis by applying the quantum 
Fourier transformation to |n). This is equivalent to the 
phase states \6) = {\/\/~2tt) E^Lo ex P (i n ®)\ n ) which are 
not orthogonal and not eigenstates of any hermitian ob- 
servable. Therefore no precise phase measurement can be 
carried out. Second, a naive extension of the sharp bin- 
ning method to the continuous case is impossible. Note 
that any coarse-grained measurement tends to lose quan- 
tum properties [14| and lead to non-tight BI tests. From 
the above results for unsharp and regional binning we 
also do not expect strong violations by these methods 
for the continuous variable system. 

Let us consider the Pegg-Barnett phase state formal- 
ism (l5| . We approximate the quantum phase by an 
orthonormal set of phase states in a s + 1-dimcnsional 
truncated space \6,k) = (1/y/s + l) En=o ex P (in6k)\n) 
where 6 k = 8 + 2nk/(s + 1) and k = 0, 1, s. Note that 
s is a cutoff parameter (assumed here to be an odd num- 
ber) and in the limit s — > 00 there exists a 9k arbitrarily 
close to any given continuous phase. The correlation op- 
erator can then be written as E(9, <fi) = n(#)<g>n(</>) using 
the phase parity operator 11(0) = Efc=o(~ 1) k)(9, k\. 

We consider a truncated TMSS \tp a ) = 

(sechr/Vl - tanh 2s+2 r) EJUo ta nh" r\n, n) which 
tends to the s + 1-dimensional maximally entangled 
state for r — » 00 and to the TMSS for an infinite cutoff, 
s — > 00. The preparation of this state can for instance 
be achieved by the optical state truncation method (l6j . 
The expectation value of the Bell operator is given by 



QM 



(iPs\Etf, 4>) + E(fi, 4') + 0) ~ E(6', ^)\^ s 



= 4V2 



tanh 2 r 
l + tanh s+ V 



(8) 



when 6 = 0, 6' = ir/(s + 1), 4> = -tt/(2s + 2), and (/>' = 
7r/(2s+2). Fig.[^a) shows its monotonic increase against 
the squeezing rate r for different cutoff parameters s. For 
any finite s and S > there exists a squeezing parameter 
r above which Bqm > 2\[2 — 5. The required squeezing 
for this violation is r > | ln[(l + f(s, 5))/(l - f(s, 6))] 

where f(s, S) = [{2y/2 - \/^V25 - 6*)/(2y/2- 5)} 2 ^ s+1 \ 



4 




tanh (r) 



FIG. 2: (a) Expectation values of the Bell operator for trun- 
cated TMSS with cutoff parameters s = 1 (solid), s = 9 
(dashed), and s = 99 (dotted), (b) The shaded region indi- 
cates the values of r for which the BI is violated; a violation 
better than £>qm > 2y/2 — 8 occurs above the curves shown 
for 8 = 0.01 (dotted), 8 = 0.001 (dashed), and 8 = 0.0001 
(solid). 



The shaded region in Fig. HJb) indicates the values of r 
for which the BI is violated £>qm > 2 and a violation 
better than Bqm > 2y/2 — 5 occurs for values of r above 
the corresponding curves for different 5. Violations ar- 
bitrarily close to the maximum value 2\/2 can thus be 
achieved by sufficiently strongly squeezed states for any 
finite value of s with r — > oo corresponding to the EPR 
state. Remarkably, this is in contrast to previous types 
of Bis which were not able to get arbitrarily close to this 
bound for the EPR state. However, we should note that 
for large s one here again faces difficulties in performing 
precise measurements due to the indistinguishability of 
two local measurements as tt/(s + 1) — > for large s. 

Finally, we discuss the relation of our optimal BI 
with the BW inequality proposed in [llj. There, lo- 
cal measurements are performed in the basis obtained 
by applying a Glauber displacement operator D(a) on 
the number states \n). The measurement basis is writ- 
ten as | a, n) = D(a)\n) with a an arbitrary complex 
number. The displaced number operator is defined as 
n a = D(a)nD* (at)- Since n a \a,n) = n\a,n), the cor- 



relation operator is given by E(a,[3) = 11(a) ® ET.(/3), 

where 11(a) = J2^Lo(~ l)"l a ' n ) ( a ; n \ ^ s the displaced 
parity operator. Using this notation the BW inequal- 
ity becomes equivalent to Eq. (J7J), which shows that it 
is a tight BI for continuous variable systems. However, 
the maximal expectation value of the BW inequality was 
shown to be 2.32 < 2y/2 [TO], while our type of BI asymp- 
totically reaches the bound 2y/2. This shows that the 
optimal measurement bases for this non-locality test are 
obtained by a quan tum Fourier transformation on the 
standard bases 17[, i.e. each of them is mutually un- 
biased to the standard basis. This may also provide a 
useful insight about the optimality of measuring in mu- 
tually unbiased bases for cases with more than two local 
measurements [l8j . 

Conclusions-We derived, for the first time, a BI in even 
d-dimensional bipartite systems which is maximally vio- 
lated by maximal entanglement and is also tight. These 
are desirable properties for Bis in high-dimensional sys- 
tems 0, Q . Our BI is found by perfectly sharp binning of 
the local measurement outcomes. It can be used for test- 
ing quantum non-locality for high dimensional systems, 
for instance it coincides with the result for heteronuclear 
molecules by Milman et al. [il^ ]. Furthermore, we ex- 
tended our studies to continuous variable systems and 
demonstrated strong violations asymptotically reaching 
the maximal bound 2y/2 for truncated TMSSs by par- 
ity measurements in the Pcgg-Barnctt phase basis. This 
provides a theoretical answer to the question of how max- 
imal violations of Bis can be demonstrated for the EPR 
states in phase space formalism Q . In the future we will 
investigate the susceptibility of violations of our Bis to 
measurement imperfections. In this context it will also be 
valuable to search for additional optimal Bis comparing 
their properties and extending optimal Bis to multipar- 
tite systems. 
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